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АВ5ТКАСТ 


Recognition of the need to develop optimum diffusers 
for advanced centrifugal compressors, resulted in the design 
and manufacture of a novel low speed test facility for 
centrifugal diffuser testing. The CDTD was designed to 
allow the flow angle and wall boundary profiles into the 
test diffuser to be controlled by variable geometry in the 
flow generator. The present study reports on the design 
of the flow generator and the analysis of the internal flow 
using a NASA computer code (MERIDL). First test results 
are given and are compared with the results of a control 
volume analysis. The flow angle control technique was 
found to work effectively but to give somewhat smaller 
angles (by 4°) than were predicted. It was concluded that 
the information obtained would allow scaling of the device, 
however an analysis code was needed which would accept the 


real physical boundary conditions. 
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I. o NERODUCITION 


Centrifugal compressors are currently used in gas 
turbine power units for turboprop aircraft, helicopters 
and auxiliary power production. The capability of achiev- 
ing high pressure ratios per stage (12:1), and the simplicity 
Gr design & fabrication, are reasons for using centrifugal 
rather than axial stages in smaller engine types. The 
main disadvantage is the progressively lower efficiencies 
obtained as the pressure ratio is increased. 

A major contributor to centrifugal compressor ineffi- 
ciency is the performance of the diffuser which closely 
follows the impeller of the compressor. The purpose of the 
centrifugal diffuser is to convert most of the kinetic 
energy of the flow entering the diffuser vanes into static 
pressure, with the highest efficiency attainable. 

The design of centrifugal diffusers is presently based 
Mainly on experimental results for two dimensional and coni- 
Gal diffusers. While numerical methods are currently under 
development, computer solutions to viscous, three dimensional, 
unsteady, transonic flows, with adverse static pressure 
gradients are not yet available as verified tools to be 
used to optimize designs. 

Centrifugal diffusers are currently tested only as 


components of high speed compressors or gas turbine engines. 
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This approach is very expensive and does not yield detailed, 
accurate information which is necessary to confirm design 
systems or to provide the basis for improved theoretical 
analysis. 

A new test facility [Ref. 1] has been proposed and built, 
the so-called CDTD (Centrifugal Diffuser Test Device), 
which has as its main purpose the large scale, controlled 
simulation of the time-averaged inlet flow to centrifugal 
compressor diffusers. The device will allow the detailed 
evaluation of proposed diffuser geometries at an acceptable 
expense, and will permit validation of new computer analy- 
sis codes for diffusers operating in a fully subsonic flow 
regime. 

This report deals only with the flow generation within 
the CDTD, the control of which was to be effected in a novel 
way. The flow was supplied through a central, swirl gener- 
ating cylinder (GC) and the average angle of the flow into 
the diffuser under test, was to be controlled by covering 
or exposing more or less of the cylinder's length. 

The design and the calculations on which the device was 
Fr G, are given in Chapter II. Іп the same section are 
also described the results obtained of modeling the flow- 
field generation using the NASA computer code MERIDL [Ref. 
21. In order to use the code, in its standard form, 
artificial boundary conditions were used to obtain nearly 


representative conditions at the physical boundaries. 
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Selected results from the initial program of measure- 
ments, Carried out on the CDTD, operating in effect without 
a diffuser, are given in Chapter III. 

The flow angle control principle, was shown to work 
qualitatively, however the angle produced for a given 
exposed length of the central cylinder, were found to be 
less (by 5°-10°) than were predicted in the design 
ealculations. 

Reasons for the observed differences between design and 
test results are discussed in Chapter IV. It is shown 
(Appendix D) that a control-volume analysis of the flow 
generation, in which the geometry of the swirl generating 
cylinder surface, mixing and wall effects are included, can 
be used to predict the measured results. 

It is concluded in Chapter V that by using the results 
of measurements to establish reasonable values for unknown 
factors in the control-volume analysis, a means of scaling 
the design is obtained. 

Finally, it is recommended that an axi-symmetric, inviscid 
analysis of the flow generation, in which the physical 
boundary conditions are more properly represented, be 


carried out. 
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шие СЕРРБЕЗЛСОЫ 


ES DESCRIPTION 

The Centrifugal Diffuser Test Device (CDTD), is located 
іп the Cascade Building (213) at the Turbopropulsion Labora- 
tory complex of the Naval Postgraduate School. A view of 
the apparatus is shown in Fig. 1. A schematic is shown in 
Ша. 2. 

The philosophy in the design of the CDTD was to provide 
a nearly tangentially-directed, uniform airflow from a large 
cylinder, have a minimum decrease in tangential velocity 
between the cylinder and the test diffuser (conserving angu- 
lar momentum) while the exposed length of the cylinder would 
control the radial component of velocity and hence the average 


flow angle (8.) entering the test diffuser. 


2 
The device consists of the following major components 
with their corresponding functions. 
езек ава Cylinder (Fig. 2, Index 5) 

The perforated cylinder serves as a major structural 
component of the CDTD, holds the north and south plates in 
parallel and concentric alignment and provides uniformly 
distributed air at low velocity around its periphery. 

EMEECcnenstesngocvwMDder (Fig. 2, Index 10) 


The generating cylinder is placed around the per- 


forated cylinder and consists of the following components. 
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Gee cieralsceetron (hig, 2; Index 12) 

The central section is made of sets of small 
nozzles (not 2-D vanes as shown in the schematic) formed 
by pressing & bending sheet metal strips and soldering them 
to form a drum. A view of the complete section is shown in 
Fig. 3. The nozzles provide an almost tangential velocity 
which is controlled by the pressure supplied at the inlet. 

Bess orehsszsautch End Cylinders (Fig. 2, Index 11) 

These support the central section. They each 
contain a short length of porous section inside of which 
cylindrical throttles (Fig. 2, Index 6) are arranged to 
slide. The throttles, operating separately, expose or cover 
the porous sections, so controlling the air supplied to the 
jet walls. 

БК е Walls (Tig. 2, Index 3) 

These are annular end walls containing nozzles dis- 
tributed radially, so that nearly tangential injection is 
produced through them (the same nozzles as for the gener- 
ating cylinder are used). By-pass air, depending on the 
throttle positioning is routed through the jet walls, to 
effect changes in the velocity distribution at the diffuser 
inlet. The jet walls are also the means for controlling 
the main airflow into the test section (TS). The mass flow 
through the swirl vanes, is controlled by moving the jet 
walls axially, and so exposing more or less of the swirl 


Vane area. Thus the velocity profile and the average flow 
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angle (85) at the test section, are controlled by the jet 
walls and associated throttle positions. 
on 29Eher Components 
The outer casing (Fig. 2, Index 9) contains a con- 
Brsetion contour (Fig. 2, Index 8), which is designed to 
turn the flow into a radial, parallel-wall test section 


containing the diffuser vanes under test. 


B. DESIGN CALCULATIONS 

The CDTD design calculations are presented in Appendix 
A. The concepts of continuity and conservation of angular 
momentum were applied to an incompressible, inviscid flow. 
The detailed geometry of the inlet surface and effects of 
friction were not considered. 

Selected flow angles (В) апа jet wall spacings (La), 
were input into the derived formulae, and the results are 


given in Table Al and Figure 29. 


C. FLOW ANALYSIS 
A FORTRAN computer program MERIDL [Ref. 2] was used to 
obtain a prediction of the flowfield generated by the CDTD. 
Specifically, in the design phase the main interest was 
to investigate the velocity behavior along the contraction 
KI Our (Fig. 2, index 8), in order to verify that the flow 
is turned without adverse deceleration ahead of the test 


vanes, into the diffuser. 
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MERIDL, besides its general applications for calculating 
flows through blade rows in turbomachines, can obtain 
solutions for flow in annular ducts without blades. 

The flow must be essentially subsonic and the solution 
obtained is for two-dimensional axi-symmetric, compressible, 
shock free flow. Upstream and downstream flow conditions 
can vary from hub to shroud, and provision is made for an 
appropriate correction for loss of stagnation pressure. 

The analysis consists of the solution of the simultan- 
eous, non-linear, finite-difference equations for the stream 
function (Y). 

The following assumptions are made: 

(1) The fluid is a perfect gas with constant specific 
heat Co 

(2) The only forces along a hub-shroud orthogonal mesh 
line are those due to momentum and pressure gradient 
(viscous forces are neglected in that direction). 

(3) There is no heat transfer. 

(4) The upstream and downstream boundaries of the 
solution region are orthogonal to the streamlines. 

(5) The stream function is zero at the hub and 1.0 at 
the shroud. 

The program generates an orthogonal mesh in the space 
between the hub and the shroud, by dividing it into equal 
increments along several hub-shroud lines. Spline curves 


are fit through the resulting points to obtain the streamwise 
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orthogonals. The normal orthogonals are obtained by a 

БЕЕСІстог-сСоттесток technique. The solution of the equations 

follows an iterative method with successive overrelaxation. 
1. Flowfield Modelling 

An equivalent flowfield to the one generated by the 
СОТО had to be formulated in order to overcome MERIDL's 
inability to accommodate the real physical boundary condition 
along the generating cylinder. 

The author of Reference 2 was consulted and the 
geometric model shown in Fig. 4 was developed after appro- 
priate calculations were carried out (Appendix B). 

The inflow stream was introduced with a constant 
whirl (R»r = const.) far upstream of the TS. The channel 
height was doubled so that the streamlines, close to the 
region of the contraction contour, approximated the actual 
Blow pattern generated by the device 

The input files to MERIDL (Appendix E) were gener- 
ated and the program was run using the Naval Postgraduate 
School's IBM 370-3033 computer. 

2. Results 

Selected tabulated results (for j = 0.0 and y = 1.0) 
from the converged output files generated by MERIDL are 
given in Appendix F. 

The TEKTRONIX 618 plotter was used to plot the 
streamlines of the computed flowfield (Fig. 5a). An 
enlargement of the area of main interest around the con- 


action contour is shown in Fig. 5b. 
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The relative velocity along the outer casing wall 
(4 = 1.0), for the selected runs, is shown plotted in Fig. 6. 
The flow appears to be accelerating in the region of the 
contraction contour (stations 14-23) and then begins to 
diffuse. 


Finally in Figures 7 and 8 the velocity and the flow 


angle variations are shown plotted for the y = 0.0 stream- 
line from the surface of the GC (к, ASS SUE E) up to- the 
TS inlet (R3 = 2.l FT). The velocity is seen to diffuse 


smoothly and the flow angle decreases to almost a constant 
value at the inlet station to the test vanes. 

A check was made of the velocity and flow angle 
uniformity at the location of the GC. The results are given 
in Table 1. Values between stations were obtained by linear 
interpolation of the tabulated output from MERIDL. The 
results show that the conditions expected, at the surface 
of uniform velocity and angle, appear to have been achieved 
with a maximum deviation of .5% in the velocity and 2.4% 
in the flow angle (81) at outlet flow angles from 60? to 
po". 

Comparisons of the predictions with measurements 
will be made when pressure taps have been installed along 
the outer casing and the contraction contour, as discussed 


M Subsection III.A.2. 
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III. PRELIMINARY TEST PROGRAM 


Preliminary measurements were carried out as shown in 
Figure 9, using a short vaneless diffuser, obtained by re- 
placing the plexiglass walls (Fig. 2, Index 13), with locally 
fabricated plywood walls (Fig. 9, Index 4). The walls ended 
at the intended location of the leading edge of the test 
vanes. 

The data were taken using two pressure probes, one United 
Sensor 5-hole and one Kiel probe. Both probes were connected 
to a manometer board, inclined at a 30° angle for greater 


sensitivity. 


A. INSTRUMENTATION 

1. Installed for Preliminary Measurements 

The instrumentation arrangement is indicated in 

Figure 9. Тһе 5-hole probe (Fig. 9, Index 1) was attached 
to the south half of the outer casing which could rotate to 
facilitate circumferential measurements. The radial position 
of the probe corresponded to the leading edge of the test 
vanes (к. = 25"). The probe could traverse axially in order 
te obtain data for the total pressure and flow angle dis- 
tributions, as well as for the velocity profile, produced 


across the outlet from the plywood walls (parallel to the 


Main axis of the device). 
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The Kiel probe (Fig. 9, Index 5) was moved radially 
from the GC surface (Ry = 19"), up to the intended location 
of the LE of the test vanes (R. = 25"), and axially (from 
wall to wall). “The purpose of measurements inside the flow 
generator was specifically to obtain: 

a. The flow mixing pattern out of the swirl vanes and 
from the GC up to the outer casing. 

p. Wake visualization close to the surface of the GC. 

c. The blockage factor (K, ) for the flow, at the exit 


Ш 
of the nozzles. 


ЖЕ Гизат тес ror Preliminary Measurements 

A set of pressure taps was designed to be drilled 
on the contraction contour and on the outer casing, as 
shown in Figures 10a and 10b in order to obtain data to 
compare with the results generated using MERIDL. 

The proposed arrangement is shown in Fig. 10a, 
drawn on the north (fixed) part of the apparatus. The 
circumferential tap distribution is indicated in Fig. 10b, 
for both rotating and fixed walls. 

A tap size of .04 inches diameter (.001 m) was 
selected. Pressure data can be recorded using the labora- 
tory's HP 3052 Data Acquisition System with Scanivalve 


Enterface. 


B. PROGRAM OF MEASUREMENTS 
The goals of the preliminary measurements were: 


a. To verify the effectiveness of the flow control method, 
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b. To visualize the flowfield, 

с. To verify circumferential symmetry, 

d. To measure the total pressure and flow angle (55) dis- 
tributions at the inlet to the test section. 

A series of tests were conducted, at a plywood wall 
Брастпа ОЕ 2.62" (varying slightly circumferentially), at 
selected jet wall spacings and throttle settings. Probe 
measurements were also taken at various circumferential 


positions. 


E RESULTS 

Data from selected runs are presented in Tables 2 through 
7. 

The data obtained from the measurements were reduced 
using the methods of averaging derived in Appendix C. 

The averages were obtained assuming incompressible flow and 
using the equations for conservation of mass and angular 
momentum. 

A computer program was developed, to carry out the re- 
quired calculations. Spanwise (gapwise) integrations were 
accomplished using the subroutine DATINT from Reference 3. 

gcom the computations, blockage factors M and mass 
averaged flow angles (B,) were obtained at the inlet of 


the test section for various configurations of the apparatus. 
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IVES DISCUSSION 


fe EXPERIMENTAL RESULTS 
The initial results showed that: 
1. Тһе flow control method was working, but the flow 
angle was not affected as much as was expected by 
jet wall spacing changes. 
2. The measured maximum average flow angle В) was less 
than the design angle of 70°. 
3. The flowfield generated by the device was qualitatively 
as expected and mixing of the jets from the swirl 
vanes occurred very quickly. 
4. Circumferential symmetry of the measured flow angle 
8, at the exit of the plywood walls was not achieved. 
In trying to explain the results, the importance of the 
leak which was forced to occur at the interface of the GC 
Surface and the jet walls, in all but the maximum (25 = 24.375"), 
spacings was examined. Table 2-9 gives data for 14 ZZ T S 
A survey with a tuft at minimum jet wall spacing Tz = 7.0") 
showed that significant axial component of velocity anda 
vortex flow were generated at the gap, as shown in Figure IV.l 
below. 
A vortical recirculation region was established near 
the jet wall, between the outer casing (higher static 


pressure region) and the generating cylinder surface (lower 
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Figure IV.l. Schematic of the Vortex Flow Generated 
ace ERE CCST гпсегтасе 


static pressure region), which would contribute to a reduc- 
tion in the total pressure and flow angle along the contrac- 
pon cone. In addition, non-uniformity in the leak 
сшесипвегепъта! Ту, would result in non-uniformity in the 
Circumferential measurements. 

Sealing of the gap, using strips of upholstery piping, 
resulted in the data given in Tables 1-17. Noted were: 

l. An increase in flow angle В) by approximately 4° 
(Figs. 11 and 12) and improvement in the total pressure 
drseribütion (Figs, 13 and 14). 

2. Satisfactory circumferential uniformity in flow angle 


(Figs. 15 and 16) and total pressure (Figs. 17 and 18). 
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3. The effect of the jet wall spacing on В) ӘНШІ 17) 
approached more closely the design expectation (Fig. 
29), the difference which remained will be discussed 
in the following subsection. 

in Fig. 20 is shown the effect of JW spacing on the 
total pressure distribution at the location where the 
measurements were taken. 

Figures 21 through 24 show the effect of the JW throttle 
setting on the flow angle (Figs. 21 and 22) and on the total 
pressure (Figs. 23 and 24), for the two JW spacings at which 
reliable (leak-free) measurements were taken 

Figures 25 through 28 show circumferential measurements 
of the flow angle and the total pressure distributions for 


the jet walls in the fully retracted position. 


BRO THEORY AND EXPERIMENT 

In Fig. 29 are shown the: results for the measured flow 
angle (averaged as described in Appendix C) and the values 
obtained in Appendix A, using the design calculation method, 
for the wall spacing at the time of the tests of 2.62". 

The experimental values are seen to be lower by 4° at maximum 
and 8° at minimum jet wall spacing. 

The design calculation was based on an idealized flow- 
field and was immediately questioned. A complete control 
volume analysis (Appendix D) was carried out in order to 
take into account the geometry of the generating cylinder 


mud the contributions of wall friction. 
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Examples of results are plotted in Fig. 29. Тһе соп- 
trol volume analysis appears to agree well with the measured 
values when the JW's are widely spaced, but to predict a 
Stronger effect of reducing the spacing than is actually 
measured. The significant differences between the idealized 
design predictions and the control volume analysis however, 
imply that the geometry of the generating cylinder surface 
and effects of friction need to be included in analyzing 
the flow. It is noted that the control volume analysis 
includes one unknown blockage factor, and one which must 
be input from measurements. Further examination needs to 
be made of the uncertainties in the analysis, in order to 
explain the differences in the slopes of the lines in Fig. 
ЕЭ. 

Finally in Fig. 29 are shown the calculated flow angles 
corresponding to the design value of the diffuser wall 
Spacing (2.0"). Clearly an increase of the wall spacing 
Results in an increase of the flow angle. Any desired 
range of flow angle can be obtained by changing the width 


(or span) of the test diffuser. 
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drawn: 


ir; 


V. CONCLUSIONS 


From the work reported, the following conclusions were 


The proposed method of flow control works 
satisfactorily. 

Sealing the jet wall gap is essential, in order to 
obtain the desired flow angles and circumferential 
flow symmetry. A method to provide sealing, at 
specific values of jet wall spacing must be designed 
in order to achieve the desired working range of 
diffuser inlet flow angle. 

A review of the control volume analysis is necessary 
in order to explain the measured flow angles. 

While results obtained with MERIDL were not unrealistic, 
the development of a computer code which will accept 
the proper boundary conditions, for the CDTD, is 
needed in order to predict more accurately the 
distributions of flow conditions at the test section 


inlet. 


An advanced code with transonic flow capability is 


required ultimately for the design of a high speed version 


er the CDTD,. 
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APPENDIX A 


CDTD DESIGN CALCULATIONS 


Assumptions: 
1. The velocity coming out of the JW is the same as at 
the GC surface (Fig. Al). 
2. No losses. 
3. Constant density (incompressible) flow. 


The velocity components are given by: 


Ур. = V; cos P; даи) 
1 
апа 
Vo. = ME sin Bi А(2) 
7 
applying continuity between stations 1 and 2 (Fig. Al). 
Vo ‚2-Ас М + Vo As = vs -Е-(2т1.) 
T 1 2 
so that 
PRs oie =! - Yo Али 
D = A(3) 
| О с 


Applying conservation of angular momentum: 
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Figure Al. CDTD Schematic for the Design Calculations 


V. R = V. R А (4) 
94 1 9. 2 
Or 
Ro 
V = У d Ааа) 
91 9 > R; 


and substituting in A(3) for Vg from Eq. A(4.1) 
1 








Substituting into Eq. A(5) the values for R Zi N, 


Ber Асу” 


апа А a linear relation is established between L; and 


JW’ 
pot 857 thus 


Li = 25.97 cot В) - 1.458 A(G) 
where data have been used for the as-manufactured test section 
span (2, = 2.61") and the short plywood walls. 

For the designed span (45 = 2"), and plexiglass walls 


installed: 
L = 19.9 cot 3, - 1.458 A(7) 


Equations A(6) and A(7) can also be solved for 8, and EEG 


thus 


oe ee А(8) 


32 1 


ВЕЕ E0251 5999073296) A (9) 


52 1 


Using Equations A(8) and A(9) the following Table Al 


is obtained: 
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TABLE Al 


FLOW ANGLES (85) AND CORRESPONDING 
JW SPACING OBTAINED FROM EQUATIONS 


65. 
p. 
Ко. 
52. 
63. 


ИЗ. 


78 


po 


A(8) AND A(9) 


* 


per 


dy (in) В) ka (ай) 
21.304 65 15.643 
15989 70 1.07 
117001 75 7.748 
29 24,375 29226 24.375 
075 18.00 62.27 16.00 
257 12.00 69.456 12.00 
. 14 8.00 74.663 8.00 
2 700 79.01 7.00 





APPEND xs 


FLOW MODELING CALCULATIONS FOR MERIDL INPUT 


Assumptions: 

1. Free vortex flow: V ES = constant (from hub to tip 
at a certain een). and V, = constant with radius 
al Constant ds 

ООО Constant density р = .076 LB /FT^. 

3. Total velocity at the inlet of the TS, V. = 200 FPS 

4. Axisymmetric flow. 

5. Flow enters from Station 1-1 (Fig. 4) with constant 


whirl. 


Flow Area at TS Inlet 


bi E 2 
A, = 2TRAL, = 1.3666 ЕТ 
VER = V, cos в, = 200 cos 8. FPS 
m = V. A = 20.772 cos В св 

2 к» : EO 
Flow Area at Inlet (Station 1-1) 
2 2 2 

Ary = т (Вс - Ва) = ЧО ЛС FT 
V, = constant 
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Applying continuity between Station l-l and at the inlet 


of the test section: 


e 


DM RUN E к; S 
so that 
> 
| ТЕВЕ 12 у В (2) 
Ы DEN Ro 


Мы К = V, R Bro) 


and solving for Vo and inserting values for R; and R^, 
1 


< 
li 


ШОЛУ, B(3.1) 
1 2 


Also, 


V = V 


' = j \ 
o > біп В) 200 Sin В) (ft/sec) B(4) 


For selected values of 8, the following Table Bl was 
constructed, using the above equations to obtain the corres- 
ponding whirl (R,V, ) and the mass flow rate (m), required 


i 
as input to MERIDL. 
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TABLE Bl 


CALCULATED INPUT PARAMETERS FOR MERIDL CODE 


В) VR, m Ve, 7952 
(Deg) (Ft/Sec) (lbs/Sec) (Ft/Sec) (Ft-/Sec) 
60 100 5225 VOD 346.4 
65 84.52 2726 191526 362.5 
70 68.4 .2206 187.94 375.9 
75 91.276 2167 193.19 386.4 
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APPENDIX C 


CALCULATION OF THE MASS AVERAGED FLOW ANGLE 


When the flow is non-uniform at the measurement plane, 
average values can be defined using the conservation equa- 
tions. For example, applying conservation of angular 
momentum between the flow entering through the swirl vanes 
and the test section inlet: 


EE V) süng.dm, - f. rjv 


2У- сіп ieu = r Va ту - rjVa m, 


where V. is the mass flow average of the tangential 


velocity component which is defined as: 


<Il 


= =]. май c(2) 


9 
m m 


Similarly conservation of mass gives: 


: 1 g 
mM = E p,V, cos 8,(2TR,)AR, я > 


в. В. "2 C(3) 


where K, іс defined as the "blockage factor" and relates 
2 _ 
the area average (Ур ) and mass average (Vp ) velocity 
2 2 


components, thus 
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К = Vp "Ув С (4) 


Mensistent with Ба. C(3) and Eq. C(4) the definition of the 


Ares averaged radial component of velocity is: 


V = =] ШЕ ар C(5) 


NJ 
© 
NJ 


The mass-averaged radial component of velocity is obtained 


вом Еа. E(3), thus 


== m 
V — а C (O) 
R, о 
Тһе mass-averaged flow angle (Bo) can be defined from: 
== “ө, 
tan 8. = — er) 
V 
2» 
so that 
V 
= EE 
B, = tan (=) С(8) 
E 
2 


шо obtain the value of B from probe survey measurements, 
and blockage factor to use in the control volume analysis 


(Appendix D), the following data reduction is required: 
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1. Calculate the area-averaged radial velocity (Va R 
2 


= 1 1 
V = | vw = = | Уу 
ES LE 


R d 2 2 


COS 8,324 C (9) 
2 


2 


where, if the flow is assumed to be incompressible 


(Pi7Pj) 5 


coc Co)‏ ج 
Ae sei tuting Eq. C(10) into C(9):‏ 


А И „5 
Ир 7 к ] (Р-Р) сов 8,47 C IT) 
2 2 2 


where 8, = 8.(2) is the measured flow angle at the exit of 


2 
the plywood walls, and (P, =P) is the difference between the 


total pressure measured by the 5-hole probe and ambient. 





2. Calculate the mass-averaged radial velocity using: 
& 
= 1 : 1 2 
о Мм = = Ј урах ei 12) 
2 m m 2 ME 0 2 
R. 2 


САЗ) 


Ди 





Calculate the mass-averaged tangential velocity using: 





V, = Е] у а = o o ао C(14) 


Calculate EE using EG. C(8). 
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APPENDIX D 


CONTROL VOLUME ANALYSIS 


l. The equation for tan > 





Figure Dl. Control Volume Analysis Schematic 


Assumptions: 


1. Constant density 01 = Po 


2. Radial inlet surface S, made up of small nozzles 


I 
distributed along the GC. 


3. The wall surface (IW) consists of: 
a. The outside area of the GC; 


b. The area of the JW faces; 


19 





с. Inside area of outer casing up to the test section 


inlet. 
EM nirorm outlet conditions (Station 2). 
5. Zero axial velocity (У, ) at Stations 1 апа 2. 


Following Reference 6, the axial moment exerted by the 
Econ the walls of the control volume (Fig. Dl) is given 
DIV ; 

M = | куу, 51п 8, dm, - | Ес въз $ dm 


DE 2 


— E -- > 
ОО Бо ӨБ PD ite) 


The force ona surface element can be written in terms 


of its components normal to and in the surface as: 


dF = - pnds - ttds D(2) 
so that 
> > > = 
J R,xdF, = Í R, x (p,nds) 
а І 


ШЕСІ continuity 
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m = f dm, = f ап D(4) 


2 
S; 5. 
so that 
| dm, = ,V, sing, к A, ZNL, D(4.1) 
S 1 
1 
and 
| dm, = P,V, cos B5 Ke (27R,) 2, D(4.2) 
S 2 
2 
EnbStrtuting Eg. D(3), D(4.1) and D(4.2) into D(1): 
M = Б.ү“ іл 6 kB. A.ZN2. + R.p.A. ZN 
DENM Г 0171971 I ел 
7 NEP 
0-8-У- віп8. cos р kB, (27R,) 2, D(5) 


The different terms in the LHS of Equation D(5), can be 
written in terms of the moments of the normal and tangential 


stress components given in Equation D(2). The only contri- 


butions to the axial moment, comes from shear stress components 


(t) on the cylindrical walls, and normal stress components 


(p ) on the non cylindrical surfaces of the swirl vanes. 
м 


Thus: 


b 


M ИК К AS жө куол 9 ds' D(6) 
aN IW а “ 
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Figure D2. Stress Components on the Swirl Vanes 


From Figure D2, 
dS = 45! біп 9 pu 


sor chat 


b 
M ЕИБ У | BRAS р(8) 
азу IW 1 zo 


Combining Equations D(5) and D(8) 


b | 
E. TRAS + ZN2, | (p„py) as = m(R,)V, sin8, - RV, sin8,) 


D(9) 


living Equation D(9) for sin Br dividing both sides by 
cos Bo, and substituting Vi and У. from Equations D(4.1) 


and D(4.1), we obtain: 
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PES See x 
Б Шр E c D(10) 


tan P = Lx کےا‎ 
ET 1 


where, the pressure term 


2824 р В 
с = — J (р-р.) -- 45 D(11) 
P mV. sinß a w d R 
1 1 1 
and the friction term 
(Teds D(12) 
mV. sinß = R 
1 d: l 


The design analysis in Appendix A is the special case 


ОГ e = С. = 0. The values of ES and C. can be obtained 


using various assumptions. 


2. Calculation of the Pressure Term 
Writing, 
b 
R x 
f ESE dS = А сув 
а R; 1 


and (P. P1) # Е(К), then 


АУ, 
C = I ча (B, 7 P1) D(13) 


p 5 : 
ту, sin 83 


where (p -p,) has a maximum value of (p, -p,), or more 
Woe 1 ty 1 
meatistically, (р, -Р.), мһеге p. is the stagnation pressure 
1 


after jet mixing. 
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In order to calculate Pez” a sudden expansion is assumed 
downstream the swirl vanes. The following formula, derived 
in Reference 4, results from mixing at constant pressure: 

Pw ~ Py ШІ el 


= = 2 (1 - 


) piu 
Амтх MIX 








where AMIX is the area of the flow when total mixing has 
occurred and the velocity is uniform. 

From experimental measurements (Fig. D3), it was esti- 
mated that total mixing occurs within .3" from the GC 
surface, so as shown in Fig. D4 the area ratio is given 
Бу: 

39916 


А1/Ауту 


hence Equation D(14) becomes 


р -р, = Fe Vy x (.4765) D(15) 


м 


EUbstitutrzng to Equation D(l3), the pressure term becomes: 


20875) 


5 D(16) 
sin Bi 


It is noted that the pressure term changes the value of 


Бап В. Бу nearly 25%. 


2 
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Total Pressure Distribution in the Flow 
Generator and Near the Surface of the 
Generating Cylinder Measured with a Kiel Probe 


Figure D3. 


s 687 — 4 
D Ж лү 
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Figure D4. A,/Ayry Area Ratio Schematic 
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Aca lculation oí the Friction Term 
The friction term is approximated using an average 


velocity Se tor the slow over the inner walls of the 


CDTD; the tangential stress is given by: 


_ d -2 
t = 5 C. О Уты D(17) 


and values are required for Me. and C... 


a. Calculation of Уту 


In terms of Kiel probe measurements, 


Py =P 
5.202) 1°? D(18) 


EE | 05937 


1 


From Fig. D3 and using Eq. D(18) with the pressure 
differences (Рк-Рд) from measurements and p - .002766 <slugs/ 
вт», the following Table Dl is obtained, for the velocity 


profile close to the surface of the GC. 


Define the blockage factor for one jet as: 


9 
| У1 (2) 42 
5 
K = .ل‎ = .. 841 D(19) 


V. = V x K = 183.45 ft/sec D(20) 
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TABLE 01 


VELOCITY PROFILE CALCULATION AT THE 
SURFACE OF THE GENERATING CYLINDER 


Point # 2(рұ-рд) (in Н 0) Vi (FPS) 
1 14.0 162.0 
2 DESO 147.7 
3 eS 147.06 
4 12:12 15125 
5 10:42 190,5 
6 24.0 212.45 
7 23 219.13 
8 2953 209733 
9 №9. 2 ве. 5 


The velocity over the outer casing walls Ve can be 
approximated from the pressure difference to atmosphere 


measured at R = 20.8" (Fig. D4), then 


2080 ا 


TU M 
en 2 (5.202)] = 176.16 FPS 


p 21) 


From Equations D(20) and D(21) the average velocity over 


the inner walls МАРИ) can be estimated: 


OW 
= . pui 20 
VW L798 FPS (22) 
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ы сашептасіоп of C, 
The hydraulic diameter [Ref. 5] is used to calculate 


the Reynolds number. The Reynolds number is given by: 
T 
V D(23) 





where 


_ 4a 
а = ج‎ D(24) 


and A denotes the cross-sectional area and C the wetted 


perimeter. 


POr the СЕТО 


A = 2%, + 42, = ДЕ) pr? 
С = 2%, БЕРЕ 5 0625 СЕТ 
hence 
d, = .324 FT D(24.1) 


For v = Du от FT“/sec, Eq. D(23) gives: 


Рс оор р(25) 
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спора! universal law of friction for turbulent 


flows through pipes [Ref. 5]: 


3 





С. - 3.59 х10. D(26) 
Eubstituting all specified terms into Eq. D(12): 
-3 A мо 
457 
Я — тиу (ти, Deze) 
sin Ва T1 Vi 
4. Calculation of Outlet Angle 
USNC EC. DIG and EG. D(27) in Ec. D(10), 
_ № Am .23875 1.795 x10 2 “iw. 2 Ати 
Бап 8 = کے‎ ×) y (1 - ——— - ——M———(——)' (291 
2 > Аті sin^8 in^8 V Am] 
1 E 1 1 
D(28) 
For the values: 
Ку = 19.1 in 
К. = 25 іп 
> > 
Ar = p. = 5.0454, an 
КВ (24 = 24.375) = .925 (average measured) 
2 
k (2, = 8.0) = .83 (average measured) 
B. 1 | 
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= _ -— 
A = И кь IA DO Ko In 


T2 2 2 2 2 
А = 218.4. + 27R (L,-2.) + 21 (в2 -RÎ) - 21 ( R2- R* ) 
IW 7 OC’ 1 2 ect | ое 
NE 
= о 
By 84 


Substituting these values into Eq. D(28), the flow angle 
as a function of the exposed length (25) of the generating 


cylinder (since Кв. is also a function of 14) is given by: 


8, = tan 11(40.78)(--2)1 D(29) 
2 51 
wom БО, 0(29), for 11 = 24.375" апа k, =“ 925: 
2 
8. (25 = 24.375) = 57.13 
A۸150, for 14 = 8.0 in and кь = OS EG. D(29) gives: 
2 
= — o 
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APPENDIX E 


MERIDL INPUT FILES 
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